arXiv: 1503.00045v2 [astro-ph.CO] 15 Sep 2015 


Draft version September 17, 2015 

Preprint typeset using DT^X style emulateapj v. 5/2/11 


PAPER-64 CONSTRAINTS ON REIONIZATION II: THE TEMPERATURE OF THE Z = 8.4 INTERGALACTIC 

MEDIUM 

Jonathan C. Pober 1,19 , Zaki S. Ali 2 , Aaron R. Parsons 2,3 , Matthew McQuinn 4 , James E. Aguirre 5 , Gianni 
Bernardi 6,7,8 , Richard F. Bradley 9,1011 , Chris L. Carilli 12,13 , Carina Cheng 2 , David R. DeBoer 3 , Matthew R. 
Dexter 3 , Steven R. Furlanetto 14 , Jasper Grobbelaar 6 , Jasper Horrell 6 , Daniel C. Jacobs 15,1 , Patricia J. 
Klima 10 , Saul A. Kohn 5 , Adrian Liu 2,16 , David H. E. MacMahon 3 , Matthys Maree 6 , Andrei Mesinger 17 , David F. 
Moore 5 , Nima Razavi-Ghods 13 , Irina I. Stefan 13 , William P. Walbrugh 6 , Andre Walker 6 , Haoxuan Zheng 18 

Draft version September 17, 2015 


ABSTRACT 

We present constraints on both the kinetic temperature of the intergalactic medium (IGM) at 
£ = 8.4, and on models for heating the IGM at high-redshift with X-ray emission from the first 
collapsed objects. These constraints are derived using a semi-analytic method to explore the new 
measurements of the 21 cm power spectrum from the Donald C. Backer Precision Array_for P robin g 


the Epoch of Reionization (PAPER), which were presented in a companion paper, Ali et al. (20151. 


Twenty-one cm power spectra with amplitudes of hundreds of mK 2 can be generically produced if 
the kinetic temperature of the IGM is significantly below the temperature of the Cosmic Microwave 
Background (CMB); as such, the new results from PAPER place lower limits on the IGM temperature 
at z = 8.4. Allowing for the unknown ionization state of the IGM, our measurements find the IGM 
temperature to be above ~ 5 K for neutral fractions between 10% and 85%, above ss 7 K for neutral 
fractions between 15% and 80%, or above k 10 K for neutral fractions between 30% and 70%. We 
also calculate the heating of the IGM that would be provided by the observed high redshift galaxy 
population, and find that for most models, these galaxies are sufficient to bring the IGM temperature 
above our lower limits. However, there are significant ranges of parameter space that could produce a 
signal ruled out by the PAPER measurements; models with a steep drop-off in the star formation rate 
density at high redshifts or with relatively low values for the X-ray to star formation rate efficiency of 
high redshift galaxies are generally disfavored. The PAPER measurements are consistent with (but do 
not constrain) a hydrogen spin temperature above the CMB temperature, a situation which we find 
to be generally predicted if galaxies fainter than the current detection limits of optical/NIR surveys 
are included in calculations of X-ray heating. 

Subject headings: reionization, dark ages, first stars — intergalactic medium - galaxies: high-redshift 


1. INTRODUCTION 

Up until very recently, observational cosmology has 
lacked a direct probe of the conditions of the Universe be- 
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tween the epoch of recombination and the birth of mod¬ 
ern galaxies. Recent observations with the Hubble Space 
Telescope have found nearly a tho usand putative galax- 
ies at a redshift of 7 or above (see|Bouwens et al.|[2014 
for a recent compilation), using the Lyman break drop- 
out technique. Narrow band searches for Lyman-a emit¬ 
ters have also proven successful at finding high-redshift 
galax ies (|Tilvi et al.|2010||Ouchi et al.|2010||Hibon et al. 
201 1|) . However, these observations still principally de- 


tect the rare, most luminous galaxies, and cannot probe 
the population of more numerous, fainter galaxies. 

One of the biggest unsolved questions in cosmology 
and galaxy formation is the issue of how early galaxies 
feed back into the Universe and influence the formation 
of the next generation of galaxies. From a combination 
of Cosmic Microwave Background (CM B; Larson et al. 
2011; Planck Collaboration et al.|2015) observations and 


quasar absorption spectra measurements at high redshift 


(Fan et al. 

2006 

McGreer et al. 2015), we know that ul 

travioiet p 

hotons trom the lirst luminous objects reion- 


ized the intergalactic medium (IGM) somewhere between 
redshifts ss 6 and 15. However, it has been determined 
that the observed high-redshift galaxy population can- 


reionization of the Universe before z = 6 ( 

Choudhury 

00 
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CM 
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CD 

Kuhlen & Faucher-Giguere 2012; Emkelstein 

et al. 2012 

Robertson et al. 2013). Therefore, to under- 
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stand cosmic reionization, heating, and other feedback 
effects on the IGM, we require a probe of global con¬ 
ditions which captures the impact of the unobservable 
low-mass galaxies. 

The 21 cm line of neutral hydrogen offers such a probe. 
By observing the 21 cm signal as a function of red- 
shift, one can potentially trace the evolution of ioniza¬ 
tion, temperature, and density fluctuations on a cosm ic 
scale (for r e views of 21 cm cosmolog y, see Furlanctto 


et al. 

Loeb 


2006} |Morales fe Wyithe]1201 0 and Pritchard & 
2012|). At present, t elescopes such as the LOw 


van 


Frequency ARray (L OFAR; Yatawatta et al. 2013 
Haarlem et al. 2013 P’ the Murchison Widefi eld Array 


(MWA; [ L onsdale et al.|2009 Tingay et al.|2013| Bowman 


et al. 2013 J 

for Probing t 


et al.|2010 


^4 and the Donald C. Backer Precision Array 

re Epoch of Reionization (PAPER; Parsons 


are conducting lengthy observational cam¬ 
paigns to detect the spatial power spectrum of 21 cm 
fluctuations from the Epoch of Reionization (EoR). 

Initial measurements from a 32-element PAPER in¬ 
strument in 2011 were recently used to place an up¬ 
per limit on the 21 cm power spectrum at re dshift 7.7 
at a wavenumber of k = 0.27 /iMpc -1 (Parsons et al. 
2014|. This upper limit was stringent enough fo l place 


constraints on the temperature of the IGM, requiring 
some mechanism for heating the intergalactic gas, and 
ruling out a universe in which the thermal evolution of 
the IGM was purely adiabatic since decoupling from the 
CMB. The goal of the present work is to expand on this 
ana ly sis by usin g the more stringent upper limit from 
Ali et al. (2015) (hereafter, “Paper I”) and by using a 
semi-analytic method to model the signal, allowing for a 
more complete exploration_of t he param e ter sp ace. This 
approach improves over |Parsons et al. (2014) and Pa¬ 
per I, which constrained the spin temperature of the 
IGM by treating the brightness contrast between the 
spin and CMB temperatures as a multiplicative scalar 
on the amplitude of an analytic “patchy” reionization 
power spectrum. We review the measurements of Paper 
I in (|2j and outline our methodology in f[3] We present 
our constraints on the IGM temperature in fj4j and dis¬ 
cuss their physical implications in (J5j We conclude in 
f|6| Unless otherwise stated, all calculations assume a flat 
ACDM universe with VL m = 0.27, = 0.73, n s = 0.96, 

cr 8 = 0.82, and h = 0.7. 


2. DATA AND MEASUREMENTS 

The power spectrum measurements presented in Paper 
I represent a substantial imp rovement over the previous 
limits from PAPER in Parsons et al. (20141. While the 


previous measurements placed an upper limit of A 2 (k) < 
(41 mK) 2 at k = 0.27 ftMpc^ 1 and 2 = 7.7, the new limits 
from Paper I are significantly lower: A 2 (fc) < (22.4 mK) 2 
over the range 0.15 < k < 0.5 /iMpc -1 at z — 8.4. This 
factor of 2 reduction (a factor of 4 in the temperature- 
squared units of the power spectrum) is the result of a 
large number of improvements over the previous anal¬ 
ysis. The measurements in Paper I come from a 64- 
element PAPER array, as opposed to the 32-element ar- 


20 http://www.lofar.org 

21 http://www.mwatelescope.org 

22 http://eor.berkeley.edu 


ray used in Parsons et al. (2014). There are also several 
significant changes to the data proce s sing a nd analysis 
compared with that in P arsons et al. (20141. First, the 
appl ication of the Om nicapqredundant calibration pack- 
to the visibilit 


age (Zheng et al. 2014) 


: visibility data substantially 
improves the calibration and reduces the variance among 
measurements from physically redundant baselines. Sec¬ 
ond, the application of optimal fringe rate filtering has 
the effect of upweighting higher signal-to-noise regions 
on the sky and can limit the contamination from fore¬ 
ground emission at the northern and southern horizons. 


optimal quadratic estimator formalism (|Tegmark 

1997 

Liu & Tegmark 

2011 

Dillon et al. 

2013) m coniunc- 

tion with the de 

ay spectrum approach (Parsons et al. 


2012b I is used to downweight spectral eigenmodes that 


show significant fore ground contam i nation . In a signifi¬ 
cant change from the Parsons et al. 2014]) analysis, only 
the covariance between frequency channels within a sin¬ 
gle baseline’s measurements is used in the weighting, as 
opposed to the covariance removal techniques that fo¬ 
cused on inter-baseline covariance. See Paper I for a de¬ 
tailed description of each of these three new techniques. 

The new power spectrum measurements from Paper I 
are shown in Figure [l] (equivalent to Figure 18 in Paper 
I). The left hand panel shows P(k) in nxK 2 (/i _1 Mpc) 3 , 
while the right hand panel plots the dimensionless power 
spectrum A 2 (k) in mK 2 . In both panels, black points 
represent the new measurements with 2<r error bars de¬ 
rived from bootstrapping, while vertical dashed black 
lines represent the nominal horizon limit, beyond which 


Parsons et al. 

2012b 1 


Pober et al. 

i the right hand paneJ 

are the ex- 


pected theoretical 2er noise limit (i.e. 95% of points 

should fall under this line if the measurements are con¬ 
sistent with thermal noise; dashed-cyan) , a m odel 21 cm 
power spectrum at 50% ionization from Lidz et al. (20081 
(magenta), and three previous upper limits on the 21 cm 
signal: the Giant M etrewave Radio Telescope measure¬ 
ment at z = 8.6 (Paciga et al. 2013[ yellow triangle); the 
MWA measurement at z = 9.5 (|Pillon et al. 2014] pur¬ 
ple triangle); and the|Parsons et ai.|(|2014|) measurement 
(green triangle). 


3. METHODOLOGY 

As a way of placing these results in the context of a 
large and uncertain parameter space, we identify two pa¬ 
rameters as being the dominant contributors to 21 cm 
power spectra with amplitudes at the level of the PA¬ 
PER constraints: the average spin temperature of the 
emitting (i.e. neutral) gas — which at these redshifts is 
set by the kinetic temperature of the gas — and the aver¬ 
age neutral fraction. To understand the physical condi¬ 
tions under which these two parameters drive the 21 cm 
power spectrum, it is worthwhile to keep the brightness 
temperature contrast between the 21 cm signal and the 


23 https://github.com/jeffzhen/omnical 
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Fig. 1.— Left: The measured P(k) in units of mK 2 (/i 1 Mpc) 3 over both positive and negative line-of-sight wavenumber fey. Right: The 

dimensionless power spectrum A 2 (k) = 2^-P(/c) in units of mK 2 versus \k\. In both panels, black points represent the new measurements 
with 2cr error bars derived from bootstrapping, while black dashed lines represent the nominal horizon limit to flat spectrum foreground 
emission. Also shown in the right hand panel are the expect ed theoretical 2 cr upper bounds of a noise-dominated measurement (dashed 
cyan), a model 21 cm power spectrum at 50% ionization from Lidz et al. ( 2008[) (mage n ta), and three previous upper limits on the 21 cm 
signal: the Giant Metrewa ve Radio Telescope meas urement at 2 : = 8.6 Paciga et al. 2013] yellow triangle); the MWA measurement at 
z = 9.5 ( [Dillon et al.|2014| purple triangle); and the |Parsons et ah] ( |2014| ) measurement (green triangle). 


CMB, (5Tb, in mind: 


21 cm power spectrum. 


ST b (u) « 9x H i(1 + <5)(1 + z ) 2 


1 - 


Tcmb(z) 

T S 


H(z)/(l + z) 
ck>||/dr|| 


mK, (1) 


where ^hi is the global neutral hydrogen fraction, z is 
the redshift, Xcmb is the temperature of the cosmic mi¬ 
crowave background, Ts is the spin temperature, H(z) is 
the Hubble parameter, and di>||/d7'|| is the gradient of the 
prope r velocity along the line of sight (Furlanetto et al. 
2006). It is worth explicitly stating that all the terms in 
Equation [l] can have different values at different spatial 
locations in the universe. When we refer to the mor¬ 
phology of, e.g., the ionization or spin temperature field, 
we are referring to the spatial distribution of the fluctua¬ 
tions in these quantities. We often quantify the statistics 
of these fluctuations in cosmological Fourier space using 
a power spectrum. If we define a fractional brightness 
temperature perturbation, 6 21 (x) = [6T b (x ) - 6T b ]/ST b , 

the power spectrum, P(fc), is given by the ensemble av¬ 
erage of the square of the spatial Fourier transform of 
this perturbation: 


(M*)M#)) = (27T) 3 S(k-k')P(k), (2) 

where the unsubscripted S is a Dirac delta function. 

Using these relations as a framework, we can now dis¬ 
cuss the impact of our two principal parameters on the 


3.1. Spin Temperature 

In the brightness temperature 5T b , the spin temper¬ 
ature enters as a ratio with the CMB temperature: 
[1 — Tcmb{z)/Ts\- If the spin temperature is much larger 
than the CMB temperature, this term saturates at a 


the spin temperature is already very large (e.g. 

hirlan- 

etto 

2006 

Pritchard & Loeb 

2008 

Pober et al. 

2014). 

It is thought that the emission of ultraviolet p 

notons 


from the first luminous objects couples the spin tem¬ 
perature to the kinetic gas t emperature field through 
the W outhuysen-Field effect (Wouthuysen 1952| Field 


1958) And, in most models of early star and galaxy 


formation, the kinetic gas temperature has been raised 
to a very high level through heating from X-rays from 
the first high-mass X-ray binaries (HMXBs). (Recent 
work by |Fialkov et al.|2014 has called this last statement 
into question, motivating the “cold” reionization scenar¬ 
ios we consider here.) However, it is clear that a very 


24 Using Equation 7 from |McQuinn & O’Learyl 
(2012), we estimate that a star formation rate density or 
2.5 X 10~ 3 MgMpc -3 yr —1 is necessary for the Wouthuysen-Field 
effect to couple the spin and color temperatures in the IGM by 
z = 8.4. The observed high redshift star formation rat e density 
is nearly an order of magnitude higher than this value (|Bouwens| 
|et al.|2014[ |McLeod et al.|2014| ), making the assumption that the 
spin temperature is equivalent to the kinetic temperature of the 
gas a valid one. 
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low value of the spin temperature — as can occur if X- 
ray heating is inefficient — can make the relevant term in 
the 21 cm brightness temperature (Equation)!]) large and 
negative, meaning the hydrogen gas is seen in absorption 
relative to the CMB. It is worth stressing that in our 
model Ts corresponds to the mass-averaged spin tem¬ 
perature. Like the ionization and density fields, the spin 
temperature also fluctuates spatially and contributes to 
the overall 21 cm power spectrum. As we discuss below, 
we do not vary the morphology of the spin temperature 
fluctuations in our simulations, but rather vary only the 
total intensity of heating, effectively scaling the global 
temperature field. 

3.2. Neutral Fraction 

The evolution of the shape of 21 cm power spectrum 
as reionization proceeds is largely driven by the evolu¬ 
tion in the neutral fraction, iehi, and its spatial fluctu¬ 
ations, Achi- Simulations have shown that neutral frac¬ 
tion largely serves as a time coordinate during reioniza¬ 
tion; put another way, the shape (and, to a lesser de¬ 
gree, amplitude) of the power spectrum can largely be 
mapped one-to-one to the global neutral fraction, inde¬ 
pendent of th e redshift at which the neutral fraction ac¬ 
tually occurs jMcQuinn et al. 2006 Lidz et al. 2008 Zahn 
et al. 2011; Pober et al. |201 -1[). (It is worth noting that 


these earlier simulations were run under the assumption 
of Tg TcmB) but we find in our current simulations 
that the neutral fraction remains the principal factor in 
determining the shape of the 21 cm power spectrum even 
when this assumption is relaxed.) 

3.3. Modelling Framework 

The basic approach of this work is to explore the 21 cm 
power spectra that occur in the two-dimensional param¬ 
eter space of (T s ,Ihi)- Paper I explores this parame¬ 
ter space using a toy, analytic model for the ionization 
fluctuation power spectrum from “patchy” reionization; 
A f{k) = (x H i - zj^)/ In (fc max /fc m in ), and scaling its am¬ 
plitude to model the effects of a. glob ally cold spin tem¬ 
perature (Parsons et al. 2012a 2014). We undertake a 
more physically motivated mapping of this sp ace, using 
the p ublicly available 21 cmF ASTp^l code vl.04 (Mesinger 
fe Furlanetto|200~7 Mesinger et aH]2011 ). Since each run' 
of this code gives a full ionization history, we need to 
vary only the spin temperature as a function of neutral 
fraction, which we accomplish by varying the X-ray pro¬ 
duction efficiency, Qx ■ By effectively lowering the num¬ 
ber of X-rays produced per stellar baryon, we reduce the 
rate at which the gas is heated, allowing us to produce 
power spectra for a large number of low spin temper¬ 
atures. Thus, we use the new PAPER observations to 
place constraints on the relative timing of the EoR and X- 
ray heating epochs: the IGM can not be too cold outside 
of the cosmic HII patches, because the resulting 21 cm 
power would be too large at z ~ 8.4. Our methodology 
strives to quantify this statement, using fiducial EoR and 
X-ray heating models (all parameters other than Qx in 
21cmFAST are kept at their default fiducial values) and 
varying their relative timing. 

In our simulations, we assume that the galactic X-ray 
luminosities follow a power law with a spectral energy 

25 http://homepage.sns.it/mesinger/DexM_21cmFAST.html/ 


index of 1.5, down to photon energies of hv > 0.2 keV. 
These values are consistent with Chan dra observa tions of 
local, star-forming galaxies (Mineo et al. 2012a b), whose 
soft X-ray luminosities (relevant for heating the IGM) 
have comparable co ntri buti ons fro m hot ISM emission 
and HMXBs. Pacucci et al. (2014) recently showed that 
the shape of the SED of early galaxies can impact the 
large-scale 21 cm signal during X-ray heating, by up to 
a factor of ~ few. As explained below however, here we 
are most sensitive to the relative timing of the EoR and 
X-ray heating epochs, which ca n impact the large-scale 
21 cm power by facto rs of ~T00 (|Christian & Loeb|2013 
Mesinger et al. 2014). 


We use 21cmFAST to calculate the ionization histories 
and associated power spectra for twelve values of Qx be¬ 
tween 0 to 2 x 10 56 (roughly 0 to 0.4 X-ray photons per 
stellar baryon). Each run produces 82 redshift outputs 
spanning the range 2 = 6.14 to z = 34.51. We then in¬ 
terpolate power spectra across the 3D space (Ts, a: hi, k). 
However, these simulations do not regularly cover this 
space; since the simulations start at such high redshift, 
there are significantly more data points at high neutral 
fractions than low. However, in the range where the 
PAPER constraints are significant, we have data points 
spaced by ~ 0.6 K in Ts and ~ 0.05 in neutral fraction. 

Slices through this space at k = 0.15, 0.25, and 
0.50 /iMpc^ 1 are shown in Figure [ 2 ] As expected, the 
brightest power spectra occur at low spin temperatures. 
The drop-off in ampltiudes at high and low neutral frac¬ 
tions is also straightforward to understand. At high neu¬ 
tral fractions, the spatial distribution of xhi is relatively 
uniform, and this term does not add appreciable power 
to the power spectrum; the predominant contribution 
comes from density (5) fluctuations, which have a smaller 
amplitude. At low neutral fractions, the lack of neutral 
hydrogen leads to an overall absence of the 21 cm signal. 
At the mid-point of reionization (xhi ~ 0.5), however, 
large ionized bubbles constrast strongly with predomi¬ 
nantly neutral regions, leading to large spatial fluctua¬ 
tions in the iehi held and a large 21 cm power spectrum 
amplitude. 

One of the motivations for using such an approach 
is the elimination of redshift as a parameter to be ex¬ 
plored)^] While one of the principal goals of early 21 cm 
experiments is to determine the evolution of the cosmic 
neutral fraction as a function of redshift, the level of 
the current PAPER upper limit does not allow for such 
an analysis. And, while the PAPER measurement is at 
only one redshift, this method allows us to use that mea¬ 
surement to place constraints on the spin temperature 
without knowing the neutral fraction at z = 8.4. The 
other main benefit of this approach is that it eliminates 
the need to run simulations which only vary what are ef¬ 
fectively “timing” parameters in the 21cmFAST code: pa¬ 
rameters that change the redshift at which specific neu¬ 
tral fractions occur, while having little to no effect on 
the shape or amplitude of the power spectrum at fixed 

26 In order to correct for the effects of redshift — 21cmFAST pro¬ 
duces ionization histories where each neutral fraction corresponds 
to a specific redshift, whereas we want to compare to measure¬ 
ments specifically at 2 - 8.4 — we scale each power spectrum 
relative to the CMB temperature at 2 = 8.4: P(k, 2 = 8.4) = 

P 2 ^sr(k,z)*( %^ 
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Fig. 2.— Power spectrum values in logio mK 2 for different combinations of Tg and .rn; at z = 8.4; the three panels show the value at 
k = 0.15, 0.25, and 0.50 /iMpc - respectively. It is clear that the power spectrum is relatively flat in k. with the brightest values occurring 
at very low spin temperatures. The sharp drop off in power above Tg = 20 K is expected; these spin temperatures are approaching the 
CMB temperature at z = 8.4, and as can be seen from Equation [T] for the brightness temperature, there is little to no 21 cm signal when 
the spin and CMB temperatures are comparable!^ 


neutral fraction. 

3.4. Effect of Other Parameters 

The exact shape and amplitude of the 21 cm power 
spectrum is the result of a rich combination of astro¬ 
physics and cosmology, which is difficult to fully map 
out in any parameter space, let alone a three-dimensional 
one. In this section, we consider the effects of parameters 
other than the spin temperature and neutral fraction on 
the 21 cm power spectrum. When the IGM is cold rel¬ 
ative to the CMB, we find that these other effects are 
sub-dominant, amounting to relatively small corrections 
to our quantitative results. The reader primarily inter¬ 
ested in the IGM temperature limits placed by PAPER 
can skip to a 

In simulations where Ts Tomb, the properties of 
the sources contributing to reionization are the dominant 
drivers of the 21 cm power spectrum and its evolution. In 
the 21cmFAST code, these parameters include £, the ion¬ 
izing efficiency of galaxies — which we find has no effect 
on the shape or amplitude of the power spectrum — and 
the minimum virial temperatures of the halos that can 
produce both ionizing or X-ray photons — which has a 
small effect on the power spectrum relative to the global 
spin temperature, but predominantly changes the tim¬ 
ing of reionization and/or heating (Mesinger & Furlan- 


neutral patches jChristian fe Loeb||2013 Mesinger et al. 
2014 Parsons et al. 2014|). We note however that the 


etto 2007 Mesinger et al.||2011 Pober et al.|2014|). Sim¬ 
ulations which spanned the range of reasonable values 
for the minimum virial temperature of ionizing halos in 


Pober et al. (2014 1 only varied the peak power spectrum 
brightness at 50% ionization between « 10 — 30 mK 2 . 
This scale is far below the power spectrum amplitudes 
achievable with a cold IGM, which can range from sev¬ 
eral hundred to several thousand mK 2 , as seen in Figure 

m 

Driving the effective dominance of the spin tempera¬ 
ture in setting the amplitude of these strong (~ thousand 
mK 2 ) 21 cm signals is the contrast of the cosmic ionized 
patches (5Tb « 0 mK) with the cold ( 5Tb > —200 mK), 

27 Some interpolation artifacts can also be seen at low and high 
neutral fractions for the larger spin temperatures; these features 
occur where the 21cmFAST simulations give few data points. These 
artifacts have no effect on the conclusions of this work, which are 
concerned with the brighter power spectra at lower Tg where the 
space is well-sampled. 


uncertainty in the precise morphological structure of the 
ionization and temperature fluctuations does quantita¬ 
tively impact our constraints. Our approach largely as¬ 
sumes that the contribution of spatial fluctuations in the 
spin temperature to the 21 cm power spectrum does not 
vary with the effectiveness of X-ray heating. One expects 
that spin temperature fluctuations would be small if 
heating is inefficient and would grow larger with stronger 
heating. We quantify the change in the shape of the 
21 cm power spectrum versus heating efficiency in Fig¬ 
ure [3] For low neutral fractions, the change is small 
(~ 10%) regardless of spin temperature. At higher neu¬ 
tral fractions, spin temperature fluctuations increase the 
21 cm power spectrum at large scales (k < 0.2 /iMpc -1 ) 
by factors of as much as 3 when heating is relatively ef¬ 
ficient and the global spin temperature is high. We also 
see a decrease in 21 cm power by factors of up to ~ 30% 
at smaller scales (k > 0.6 /iMpc -1 ). However, in the 
regime constrained by the PAPER measurements (scales 
0.2 hMpc^ 1 < k < 0.5 hMpc -1 , and, as shown subse¬ 
quently, spin temperatures less than 10 K, demarcated by 
the dotted lines), we find spin temperature fluctuations 
introduce typical changes of ~ 10% or less. Therefore, 
our approximation of a spin temperature field morphol¬ 
ogy independent of heating efficiency should introduce 
only small uncertainties into our constraints. This result 
largely confirms the findings of [Pritchar d fe Furlanetto] 
(2007), and validates the approach of Paper I to use the 
global spin temperature as a multiplicative scalar for the 
overall power spectrum amplitude. 

Our methodology also implicitly assumes that the con¬ 
tribution to the overall power spectrum from the den¬ 
sity, velocity, and temperature fluctuation terms that 
contribute to A 2 (k) changes minimally over the redshift 
range spanned by the simulations. The ionization his¬ 
tory simulated by 21cmFAST is in principle independent 
of the heating history; in our simulations, the IGM is 
90% neutral at z = 12.5, 50% neutral at z = 9.5, and 
10% neutral at z = 8. In extrapolating all neutral frac¬ 
tions to z = 8.4, we have assumed that the ionization 
field provides the dominant contribution to the power 
spectrum across this redshift range, and that the frac- 
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k[hMpc x ] k[hMpc *] fc[/iMpc x ] 


Fig. 3. — The fractional change in the 21 cm power shape caused by spin temperature fluctuations as a function of global spin temperature 
(i.e. heating efficiency). The three panels show different neutral fractions xjji = 0.25 (left), xm = 0.5 (center), and :r h i = 0.75 (right). In 
the regime where the PAPER constraints apply (0.2 /iMpc -1 < k < 0.5 hMpc 1 . Tg < 10 K, marked with the dotted lines), varying spin 
temperature fluctuation morphology changes the power spectrum shape by ~ ten percent or less. 


tional contribution of the density and velocity fields to 
the power spectrum evolves relatively slowly. This is 
in general a good assumption for neutral fractions be¬ 
tween ss 10% and 90%, which are achieved over a nar¬ 
row redshift range, but makes the interpretation of high 
and low neutral fractions (where our constraints are the 
poorest) more questionable. We do expect the spin tem¬ 
perature fluctuations to grow more important at lower 
redshifts (the temperature of an overdense region grows 
faster than its density), but this is also a relatively small 
effect over the redshift range we extrapolate from. 

One additional free parameter in 21cmFAST is the mean 
free path of ionizing photons through the IGM, which pri¬ 
marily accounts for the unresolved, self-shielded pockets 
of neutral gas that limit the extent of HII regions. This 
parameter has be en sho wn to alter the shape of the 21 cm 


power spectrum 

Sobacchi & Mesinger 2014 

Pober et al. 

2014 Greig & M 

lesinger 

2015|), but principally only on 


the largest scales, at which the PAPER measurements 
are limited by residual foreground emission. 


These caveats do suggest that our quantitative results 
should not be too strictly interpreted, as we have ne¬ 
glected several effects that could change the constraints 
by ~ tens of percent. Given the scale of the current 
PAPER upper limit and the range of k modes mea¬ 
sured, however, working in the two-dimensional param¬ 
eter space of spin temperature and neutral fraction re¬ 
mains a well-motivated approach. 


4. RESULTS 

At each position in the (Ts,Xhi) space plotted in Fig¬ 
ure [2j we calculate the probability of getting the mea¬ 
surements shown in Figure [I] given our model 21cmFAST 
power spectrum at those values of (Ts,iehi)- We calcu¬ 
late the joint likelihood across all values of k measured 
by PAPER. As described in Paper I, the 2 a error bars 
plotted in Figure [l] are calculated from bootstrapping; 
here, we assume they follow a Gaussian distribution to 


allow for analytic calculation of the likelihood. We also 
make the conservative choice to treat all our measure¬ 
ments as upper limits on the 21 cm signal so that we 
only exclude models which predict more power than we 
observe. The PAPER measurements clearly detect non¬ 
zero power at several wavemodes; if treated as detections 
of the 21 cm signal, these points would exclude, e.g., the 
model power spectrum in Figure [l] for being too faint. 
Therefore, when calculating the likelihood that our data 
are consistent with a given model, we exclude any con¬ 
straints from points brighter than the model prediction. 

The constraints produced by this analysis are shown 
in Figure 0J As expected, our measurements are incon¬ 
sistent with very low spin temperatures, as these mod¬ 
els produce the brightest power spectra. The exact spin 
temperatures ruled out by our data depend somewhat on 
the (currently unknown) neutral fraction in the IGM at 
2 = 8.4. For neutral fractions between 10% and 85%, 
we can rule out spin temperatures below sa 5 K at 95% 
confidence. By narrowing the range of neutral fractions, 
our constraints grow more stringent: if the Universe is 
between 15% and 80% neutral a,t z = 8.4, we rule out 
spin temperatures below ~ 7 K, and for neutral frac¬ 
tions between 30% and 70%, we require Tg to be greater 
than « 10 K. The explanation for the poorer constraints 
at the highest and lowest neutral fractions is straight¬ 
forward: spatial fluctuations in the ionization field are 
small at the beginning and end of reionization, lower¬ 
ing the amplitude of the power spectrum, and allowing 
for a colder IGM to still be consistent with the data. If 
there is no heat injection into the Universe, cosmological 
adiabatic cooling brings the gas temperature to 1.18 K 
at z = 8.4 (if thermal decoupling of the gas occurs at 
at z = 200). Assuming the Wouthuysen-Field effect has 
efficiently coupled the spin temperature of the hydrogen 
to the kinetic temperature of the gas, our measurements 
require a gas temperature ~ 5 — 10 times larger than the 
minimum allowed by adiabatic cooling. 

Comparing these results with those from the ana- 
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Fig. 4.— Constraints on the IGM spin temperature as a function 
of neutral fraction based on the 2<r upper limits from the PAPER 
measurements; regions excluded at greater than 95% confidence 
are shaded in gray. Plotted is a slice through our 3D (T s ,x Hii&) 
space at the k = 0.25 /iMpc“ 1 . but the constraints are calculated 
from the joint likelihood across all k modes measured by PAPER. 



lytic toy model calculations from Paper I, we see gen¬ 
eral agreement, although our constraints are somewhat 
stronger. The calculation in Paper I does not include any 
of the physical effects of reionization; rather, it assumes a 
flat power spectrum between a minimum and maximum 
wavenumber and uses the integral constraint for a patchy 
reionization that / dlog/c A f{k) = Xm — x ^j. The min¬ 
imum and maximum wavenumbers roughly correspond 
to a range of ionized bubble scales; as they grow closer 
together, more power accumulates within a narrow range 
of wavenumbers. While it is difficult to directly compare 
this simple model to our semi-analytic results, it is re¬ 
assuring to find results constraining a similar range of 

Ts- 


5. DISCUSSION 

It is interesting to compare our constraints on the tem¬ 
perature of the z = 8.4 IGM with theoretical predictions. 
As stated previously, most models of early galaxy forma¬ 
tion and reionization predict fairly efficient heating, such 
that the temperature of th e IGM is much greater than 
that of the CM B by z < 10 (Furlanetto 2006 McQuinn & 
Q’Leary| |2012). Here, we consider two models that have 
the potential to result in low amounts of cosmic heating: 
an observationally based model, where we consider the 
heating produced by the currently observed high redshift 
galaxy population, and a more physical model, where 
we compare with predicti ons using the recently pr oposed 
reionization model from Robertson et al. (2015 1 (h ere¬ 
after , R15). These two models are discussed in 1 5.2 and 
§5.3[ respectively. 

We can estimate the heating we would expect by solv¬ 
ing the differential equation demanded by energy conser¬ 
vation in an expanding universe: 


d T k 
d t 


-2 H(z)T k + 


2 6 A' 

3 ksn ’ 


(3) 


where Tk is the kinetic temperature of the IGM, hs is 
Boltzmann’s constant, n is the number density of neutral 
gas particles, and ex is the energy injected in to the IGM 
per second per unit volume by X-ray sources (Furlanetto 


et al. 2006). To estimate ex, we use the star formation 
rate densities (psfr) measur ed from the observed high 
redshift galaxy popul ation by Bouwens et al. (20141 and 


McLeod et al. 
star formation 



and the local correlation between 
■ray luminosity: 


L *= 3 - 4 xl 04 °^I^ergs 


-1 


(4) 


where SFR is the star formation rate, and / y is an un¬ 
known high redshift normalization factor ( Grimm et al. 
2003 Ranalli et al.|2003 Gilfanov et al.|2004| ). Equation 
lean be related to e* by e x = 3.4 x 10 4U /x/absPSFR, 
where / a b s is the fraction of total X-ray emission that de¬ 
posits heat into the IGM. Given fx, /abs, and psfr(~), 
we can solve Equation [3] to predict the temperature of 
the IGM. 


5.1. Parameter Uncertainties 

Each of the three parameters described above is not 
well-determined at high redshift. Here we discuss the 
uncertainties in each before constructing models to span 
the range of uncertainties. 


• fx ■ Early measurements of the local star forma¬ 
tion rate/X-ray luminosity correlati on yielded, a 
value of 3.4 x 10 4o /a~ nyf™' 1 er § S_1 (Grimm et al. 


2003} GilfanoveFaL]|2004]); fx is a correction fac¬ 


tor relative to these initial measurements. More 
recent measurements ha ve suggested that the local 
fx va lue is closer to 0.2 (Lehmer et al.|2010 Mineo 


et al.pOlT 2012a|b ), while t he high-redshirt value 


still has large unce rtainties (Dijkstra et al. 2012). 


Fragos et al. (2013) model the redshift evolution of 
the H1V1XB population, and find that fx may be up 
to 10 times higher th an loc al, i. e., fx = 2.0, while 
observations from Cowie et al. (2012) find little to 
no redshift evolution in f \ up to 2 ss 6. In this 
work, we consider fx values which span the range 
0 . 2 - 2 . 0 . 


• /abs- X-rays heat the IGM by first photoioniz- 
ing a neutral atom. The resultant fast electron 
then deposits a fraction of its energy into heating 
the IGM gas itself. The relative fraction of this 
ionization energy which goes into heating depends 


tion 

of the hydrogen gas (Shull & van Steenberg 

(1985 

Chen & Kamionkowski 2004; Valdes & Fer- 

rara 

2008} Furlanetto & Stoever 2010). However, 


of the predominantly neutral regions in the IGM 
(as this gas is responsible for the 21 cm signal), the 
global ionization fraction is not the relevant quan¬ 
tity for determing the fraction of energy that goes 
to heating the gas. Rather, it is the small ionization 
fraction produced by X-rays which have penetrated 
into the neutral regions; in the limit of very small 
ionization fractions, a constant value of / a b s = 0.2 



























































































Fig. 5. — The dependence of evolution of the IGM temperature as a function of redshift on the uncertain parameters in the calculation: 
fx (left), fabs (middle), and a (right). In all panels, the solid pink curve plots our fiducial model with fx = 0.2, / a b s = 0.2, and a = —3.6. 
The dashed vertical red and blue lines show the redshift of the PAPER measurement (z = 8.4) and the Planck maximum likelihood redshift 
for an instantaneous reionization (z = 8 . 8 ^J - 4 ), respectively. The black dashed curve shows the minimum gas temperature from adiabatic 
cooling, and the black dot-dashed curve shows the CMB temperature. All predictions are based on the assumption that the observed z = 7 
galaxy population is providing the X-rays that heat the IGM. 


is a good approximation (Valdes & Ferrara 2008) 

Recently, Fialkov et al.|(|2014|) suggested that previ¬ 
ous heating calculations assumed too soft a spectra 
for X-ray emitters (specifically, HMXBs) at high 
redshift. Using a spectral model with fewer soft 
photons, they find the absorption of X-rays by the 
IG M can be reduced b y as much as a factor of 
5. |Pacucci et al.| ([2014]) argue the opposite: that 
the early galaxies have softer X-ray spectra due 
to: (i) a contribution from the thermal emission 
from the hot interstellar medium (ISM), which is 
locally found to be comparable at the relevant en¬ 
ergy range to that of the HMXBs (only the latter 
is used to motivate the scaling in Equation El); and 
(ii) lower column densities and metallicities of early 
galaxies compared to local ones, resulting in more 
soft X-rays escaping into the IGM. For robustness, 
we explore a broad parameter space, considering 
values of / a b s spanning the range 0.04 — 0.2 in our 
calculations. 


• Star Formation Rate Density. There has been some 
debate in the literature as to the evolution of the 
star formation rate density at high redshift. Mod- 



the best fits come from a more complex functional 
form with a slightly steeper slope). We consider 
models using both these power law indices, al¬ 
though for our fiducial model we choose a = —3.6, 

28 The fact that we are considering cold reionization scenar¬ 
ios ensures the amount of ionization caused by X-rays must be 
small. For predominantly neutral gas, the ratio of X-ray energy 
injection contributing to ioniza tion is no more than ~ 3 times 
that contributing to heating ( [Valdes fe F errara 2008). An in¬ 
jection of ~ 10 -2 eV/atom as heat would raise the gas tem¬ 
perature to ~ 100 K, well outside of the cold IGM regime con¬ 
strained bv the PAPER measurements. Therefore, a maximum of 
~ 3x 10 —2 eV/atom could have gone into ionizations, which could 
produce an ionization fraction no larger than ~ 2 x 10 —3 . 


since a larger high redshift psfr seems necessary 


to prod uce t he Planck optical depth (Planck Col¬ 


laboration et al. 2015 R15). To set the overall 


scale of this power law, we set the star formation 
rate density at z = 7 to be 10 ~ 2 M ( 7 )M pc~ 3 yr ~ 1 , 
consis t ent w ith t he observed values from IBouwensl 
et al. ( 2014| ) and McLeod et al. (2014) (although 


a correction for extinction does result in slightly 
higher values). 


The effect on the heating history from varying each pa¬ 
rameter is plotted in Figure [5] the left hand panel shows 
the effect of changing f x , the middle / a b s , and the right 
a. The parameters not being varied are held at fiducial 
values of fx = 0.2, / a b s = 0.2, and a = —3.6; in each 
panel, this fiducial model is plotted in solid pink. Also 
plotted is the expected gas temperature in the absence 
of any heating (black, dashed) and the CMB tempera¬ 
ture (black, dot-dashed). Dashed vertical lines show the 
redshifts of the PAPER measurement (z = 8.4; red) and 
the Planck maximum likelihood instantaneous reioniza¬ 
tion redshift (z = 8 . 8 ^ 4 ; blue). The CMB temperature 
at z = 8.4 is 25.4 K; if the gas temperature is below the 
CMB temperature, then the 21 cm signal will be seen in 
absorption. We see that relative to our fiducial model, 
only an increased value of f x heats the gas well above 
the CMB temperature, while choosing the minimum val¬ 
ues for either / a b s or a can result in an IGM temperature 
well below the « 10 K lower limit at z = 8.4 from the 
PAPER measurements. 


5.2. Comparison with Observations 

To better explore the parameter space, we plot the 
predicted z = 8.4 IGM temperature produced by the ob¬ 
served galaxy population versus both a and f x in Figure 
[6j For the fiducial value of / a b s = 0.2 (left), we rule out 
very little parameter space, but do disfavor the combi¬ 
nation of a steep a and low f x . In general, however, the 
observed galaxies heat the IGM beyond both the realm 
ruled out by the PAPER measurements and the CMB 
temperature of 25.4 K. With a lower value of / a b s = 0.04, 
chosen to represent the harder X-ray spectra suggested 
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Fig. 6 .— Predicted IGM temperature produced by the observed galaxy population as a function of the model parameters a and fx- The 
hashed cyan regions are excluded based on the PAPER constraints that T \q m > 10 K (assuming the neutral fraction is between 30% and 
70% at 2 = 8.4). Left: / a | )S = 0.20. With the fiducial absorption coefficient, very little parameter space is ruled out. However, we disfavor 
a combinatio n of a steep slope to the PsfrG) relation and lower value of fx close to the locally measured value. Right: /) 1 | JS = 0.04, as 
suggested by|Fialkov et al. |2014|). A considerable amount of parameter space in the observed galaxy model is ruled out, disfavoring the 
local value of fx = 0.2 and the steeper slope of psfr(U of a = —10.9. Even with this very low value of / a b s , there is still room for the 
observed galaxies alone to heat the IGM well above Tqmb- 


by Fialkov et al. (2014), we rule out a significant region 
of parameter space (right panel of Figure [6j) In this 
model, both the steepest slopes for models of Psfr and 
the smallest values of fx are excluded. However, there 
is still a large range of parameter space where the IGM 
temperature is heated above the limits set by PAPER. 

Of course, the constraints we have placed on the pa¬ 
rameters in our heating model can still be avoided by 
invoking other sources of heating — whether from an 
additional population of fainter, yet-undetected galaxies, 
or from other mechanisms such as active galactic nuclei 
and/or quasars ( VolonterifcJ Qncdin 2009) , shock heat ing 


(Gnedin & Shaver 2004| McGuinn & O’Leary 2012|) or 
even dark matter annihilation fEvoli et al.|2U14 ). It is of- 


necessary to reionize the IGM ( 

Choudhury et al. 

2008 

Kuhlen & Faucher-Giguere 

2012 

Finkelstem et al. 

2012 

Robertson et al.||2013| R15, 

; our measurements also im- 


IGM if X-ray heating turns out to be on the inefficient 
side of the currently allowed parameter space. The fact, 
however, that an additional galaxy population is required 
to complete reionization within the observationally al¬ 
lowed redshift range makes the constraints arising from 
the observed galaxy heating model used in this section 
relatively unphysical. If there are no fainter galaxies than 
those observed, it is implausible for the Universe to be 
significantly ionized at z = 8.4 (although high values 
of the ultraviolet photon escape fraction from galaxies 
at high redshift could allow for reionization to be com¬ 
pleted by the observed galaxy population). And, if the 
IGM is significantly neutral, the PAPER lower limit on 
Tg becomes much weaker (c.f. Figure |4|. However, if 
we remain agnostic about the source of ionizing photons, 
the PAPER measurements do place constraints on the 
sources heating the IGM, ruling out a range of models of 
X-ray heating. 

5.3. Extrapolating the Luminosity Function 


As noted above, galaxies fainter than those currently 
detected at high redshift are expected in most theoretical 
models, and are in fact necessary to complete the reion- 


contributors of ionizing photons; |Choudhury et al. 

2008 

Kuhlen & Faucher-Giguere) 

2012[ Finkeistein et al. 

2012 

Robertson et al. 2013 R15 

. While the analysis ii 

i i 5.2 


by only the observed galaxies, it is useful to explore a 
more self-consistent model, where we include the heat¬ 
ing that would be generated by the galaxies necessary to 
reionize the Universe. As a model, we use the maximum- 
likelihood star formation fit from R15, which uses the 
f our-p arameter fitting function of [Madau fe Dickinson| 


(20141: 


Psfr(-) = a p 


(1 + z) b 


1 + [(1 + z)/cp] d p ' 


(5) 


with best fit values of a p = 0.01376 ± 

0.001 MQyr _1 Mpc _3 ^fep = 3.26 ± 0.21, c p = 2.59 ± 0.14, 
and d p = 5.68 ± 0.19 Fa These valu es fit the compilatio n 
of sta r formation rate densities in IMadau fc Dickinsonl 
(20141, but where each measured Psfr value has been 


increased to include an extrapolation to fainter galaxies. 
Specifically, R15 take the measured galaxy luminosity 
function at each redshift and integrate down to 0.001L* 
to produce a total estimate of star formation rate density 
at that redshift (where L* is the characteristic luminosity 
of a galaxy). The value of psfr at z = 7 predicted by the 
R15 model is not significantly higher than the currently 
observed rates: 0.0201 q;qq 25 M Q yr -1 Mpc -3 (compared 
with the value of 0.01 Mgyr _1 Mpc _3 currently observed 
near this redshift), suggesting that half the galaxies 
necessary for reionization have already been observed. 
With fainter galaxies included in calculations of the IGM 
ionization history, the R15 model completes reionization 

29 Note, however, that there exist correlations between these er¬ 
rors, and as such, the range of star formation rate histories allowed 
by R15 is smaller than one might naively calculate. 
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Fig. 7.— Predicted heating of the IGM from the R15 best fit 
star formation rate density evolution; our fiducial model is plotted 
in solid pink. Shaded regions around each curve correspond to 
1<t uncertainties calculated from the 68 % confidence regions in the 
R15 model. Even when the parameters are tuned for the weakest 
heating of the IGM, the model predicts gas temperatures above 
those ruled out by the PAPER lower limit. 


by z sa 6, and produces an optical depth to rei oni zation 


consistent with the Planck value (Planck Collaboration 


et al. 20151. The reionization history predicted by this 


model reaches 50% ionization at z ~ 7.5; at z = 8.4, 
the redshift of the PAPER measurements, this model 
predicts the Universe is ~ 70% neutral. See Figure 3 of 
R15 for the full predicted ionization history. 

Figure [7] shows the IGM heating histories calculated 
from the R15 star formation rate density evolution model 
for all combinations of fx and / abs (a is not a free pa¬ 
rameter of the model.) We see that even when tuned for 
the weakest heating (f x = 0.2, / abs = 0.04), the R15 
model produces enough X-ray photons to bring the gas 
temperature above t he lower limit of 10 K set by PA¬ 
PER. (As argued in Fialkov et al. (2014), this model 
does still produce a relatively cold reionization, with 
Tgas ~ 7c mb at z ~ 8.) Our fiducial heating model 
of f x = 0.2, / abs = 0.2 heats the gas to ~ 80 K at 
z = 8.4. While this is perhaps somewhat cooler than the 
Ts Tcmb regime, this contrast still amounts to only 
a ~ 30% decrease in the 21 cm brightness temperature, 
and will be even less of an effect at z « 7 where the R15 
model predicts 50% ionization. 


6 . CONCLUSIONS 

We have interpreted new 21 cm power spectrum mea¬ 
surements from PAPER with a semi-analytic modeling 
framework. Using 21cmFAST to cover the parameter 
space of cold reionization scenarios (i.e. where the hy¬ 
drogen spin temperature is significantly below the CMB 
temperature during reionization), we find that power 
spectra with amplitudes above ~ 100 — 1000 mK 2 are 
generically produced for a wide range of neutral fractions, 
effectively independent of the other physical parameters 
during reionization. We cover the 2D (Ts, £hi) space 
with a suite of simulations, and find that the PAPER 
measurements rule out spin temperatures below w 5 K 
for neutral fractions between 10% and 85%. More strin¬ 
gently, if the Universe is between 15% and 80% neutral at 
z = 8.4, we rule out spin temperatures below « 7K, and 


for neutral fractions between 30% and 70%, we require 
Ts to be greater than ss 10 K. Given the recent mea- 



ization, R15 predict the Universe to be ss 70% neutral at 
z = 8.4, which would slightly lower our bounds on Ts-) 
We also explore a range of models for predicting the 
amount of heating the observed high-redshift galaxy pop¬ 
ulation provides by z = 8.4. We find that the observed 
galaxy population can generally heat the gas above 10 K 
and, thus, is not constrained by the PAPER measure¬ 
ments. However if the star formation rate density psfr 
drops off very steepl y w ith redshift (as su ggested by e.g. 


Bouwens et al. 2014 and Oesch et al. 2014) and the corre- 


lation between X-ray luminosity and star formation rate 
is close to the locally observed value of fx = 0.2, we 
find that additional heating of the IGM beyond that 
provided by the observed galaxies is required. If the X- 
ray emission from high-redshift galaxies has fewer low- 
energy photons than expected and the IGM is heated 
less-efficiently, our constraints tighten considerably, rul¬ 
ing out all scenarios with either a steep redshift evolution 
i n Psfr or a low value of fx- 
Lastly, we considered the predicted X-ray heating of 
the IGM that would be expected under the reionization 
model of Robertson et al. (2015). This model produces 
an optical depth to the CM 1 1 m good agreement with the 
new Planck measurements (Planck Collaboration et ah 
20151, and completes reionization by z « 6. We find that 


tor all combinations of heating parameters, this model 
heats the IGM above the lower limit set by PAPER and, 
for our fiducial heating model, results in a spin temper¬ 
ature well above the CMB temperature. Therefore, the 
current lower limit on Ts set by PAPER is consistent 
with the R15 model, but does not require more high red¬ 
shift star formation than that suggested by galaxy and 
CMB observations alone. 

The potential for future 21 cm studies is also clear from 
Figure [4j With each increase in sensitivity, each new 
21 cm measurement will rule out more of the (Tg,aiHi) 
parameter space. Figure [7] suggests that a spin temper¬ 
ature lower limit of ss 15 K would begin to rule out low 
efficiency X-ray heating models otherwise consistent with 
the R15 star formation rate density history; comparison 
with Figure [4] shows that placing a lower limit of this 
scale should require our a power spectrum spectrum up¬ 
per limit to improve by only a factor of ~ 2 (in units of 
mK 2 ). And, with an order of magnitude more sensitivity 
(again, in power units — only a factor of ~ 3 in bright¬ 
ness sensitivity), 21 cm studies will begin to constrain the 
properties of the first galaxies (|Pober et al.| 2014 Greig 
&; Mesinger 20151. PAPER has collected data with its 
full 128-element array for two seasons (double both the 
amount of data and number of antennas used to produce 
the constraints in this work), a nd the Hy dro gen E poch 
of Reionization Array (HERA; Pober et al. 2014) has 
broken ground and will begin operations in a lew years. 
With the continuing increase in the sensitivity of 21 cm 
measurements, we can expect to learn much more about 
the high-redshift universe in the near future. 
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